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ABSTRACT 

\ We present numerical simulations for a possible synthesis imaging mode of the Space 

Q\ • Interferometer Mission (SIM). We summarize the general techniques that SIM offers to perform 

imaging of high surface brightness sources, and discuss their strengths and weaknesses. We 
describe an interactive software package that is used to provide realistic, photometrically correct 
estimates of SIM performance for various classes of astronomical objects. In particular, we 
simulate the cases of gaseous disks around black holes in the nuclei of galaxies, and zodiacal 
dust disks around young stellar objects. Regarding the first, we show that a Keplerian velocity 
gradient of the line-emitting gaseous disk — and thus the mass of the putative black hole — 
, can be determined with SIM to unprecedented accuracy in about 5 hours of integration time 

for objects with Ha surface brigthness comparable to the prototype M 87. Detections and 
observations of exo-zodiacal dust disks depend critically on the disk properties and the nulling 
^-j- \ capabilities of SIM. Systems with similar disk size and at least one tenth of the dust content of 

C^) ' (3 Pic can be detected by SIM at distances between 100 pc and a few kpc, if a nulling efficiency 

of 10~ 4 is achieved. Possible inner clear regions indicative of the presence of massive planets 
\ can also be detected and imaged. On the other hand, exo-zodiacal disks with properties more 

similar to the solar system will not be found in reasonable integration times with SIM. 

Subject headings: space vehicles — techniques :interferometric — methodsmumerical - 
' circumstellar matter — galaxies: active 

> ' 

k> 1. Introduction 

la: 

For certain areas of current astronomical research, such as the detection and mass determination of 
black holes in galactic nuclei, spatial resolution is the main limitation. For others, e.g. the case of extrasolar 
planets, the challenge is to detect very faint signals in the presence of an overwhelmingly bright star. 
Interferometry from space has the potential to overcome both of these obstacles: it has no fundamental 
limit to spatial resolution, and interferometric nulling is by far the most powerful method to suppress 
unwanted light from a bright point source. 

The Space Interferometer Mission (SIM) will be the first mission to demonstrate the feasibility of 
this new and exciting technology. Scheduled for launch in 2005, SIM is the first new mission of NASA's 
ORIGINS program after SIRTF. It is designed primarily for astrometry at the /iarcsec-level. However, 
SIM may also have significant capabilities for synthesis imaging. In terms of spatial resolution, SIM with a 
maximum baseline of 12 m promises at least a factor of six improvement over the Hubble Space Telescope 
(HST). 
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In this paper, we discuss the potential of SIM to significantly improve upon current observational limits 
in two specific astronomical cases, namely line-emitting gas disks rotating around black holes in galactic 
nuclei, and continuum zodiacal dust disks around nearby stars. In the next section, we will describe how 
SIM performs synthesis imaging, and discuss some aspects and limitations that will affect the observations. 
Section || gives an overview of the capabilities of the code which we use to simulate various classes of model 
sources. Input models and first results of simulations for the two examples mentioned above are described 
in Sections |i| and [s]. We discuss and summarize our findings in Section |^. Appendices describe the basics of 
intcrferometric imaging, and the simulation of measurement errors in our models, respectively. 



2. Synthesis imaging with SIM 

As of November 1998, certain aspects of SIM design which are important for synthesis imaging have 
not yet been finalized. The most important of these is to have a sufficiently complete range of of baseline 
spacings available. For the purpose of this paper, we assume that SIM will have baselines between 0.5 



and 12 m in increments of 0.5 m, similar to the original design (see Allen et at 1997 for an overview). 
With these capabilities, SIM will be well suited to perform reliable synthesis imaging for a wide range of 
complicated sources. 

The basic "measurement" for synthesis imaging consists of determining the amplitude and phase of 
the interference fringes for a given baseline and spacecraft orientation. These two quantities together form 
the "complex visibility". In short, the signal of the two interfering beams on the detector depends on the 
path difference or "delay" between them. The interference pattern is a cosine wave whose amplitude and 
phase are the desired observables. If the path delay is accurately known and controllable, the complex 
visibilty can be determined by "scanning" the delay. Repeating this measurement for a large number of 
baselines then allows a more or less complete reconstruction of the sky brightness distribution, depending 
on its complexity. In Appendix A we describe in more detail the operation of the interferometer in terms of 
Fourier optics. Here, we concentrate on the performance of SIM for astronomical targets. 



2.1. Observing procedure 

1. The relative length of the baseline vector is monitored and controlled with an internal laser metrology 
system, while its exact length and orientation are determined by observing nearby bright stars. 
Pointing stability is achieved with two additional guide interferometers. 

2. The science interferometer is aligned such that the phase center (zero delay, or white light fringe) 
coincides with the position of a reference star. 

3. The delay lines are adjusted such that the phase center moves to the science target which has to lie in 
the 15° field of regard of SIM. In the case of astrometric measurements of a point source position, the 
only interesting observable is the delay needed to maximize the white light fringe. This determines 
the position of the source with respect to the reference star. Identification of the white light fringe is 
unambiguous since the light is spectrally dispersed: only the central fringe will show no phase shift 
for all spectral channels. 



4. 



In the case of an extended source that is to be imaged, the complex visibility is measured with 
the delay modulation method ( Shao fc Staelin (1977) ) or some variation thereof, which we describe 
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in Appendix B. Since the signal is spectrally dispersed, the complex visibilities are measured 
simultaneously for all spectral channels. This is the basis for the spectral synthesis technique that we 



describe in Section 2.4. 



Baseline length, orientation, and wavelength determine the (u, ^-coordinate to which the complex 
visibility is assigned. The wavenumbers u and v are defined as u = B x /X and v = B y /X, where B x 
and By are the baseline projections onto the x and y axes of the reference frame. The measurement is 
repeated for as many baselines as necessary to sufficiently fill the (w,w)-plane.0 The image can then 
be recovered from the [u,v )-data by Fourier-inversion, and spurious responses removed by computer 
restoration. 



2.2. Instrumental Parameters 

2.2.1. Spatial resolution 

The spatial resolution of a synthesis imaging observation is not a priori defined. A diffraction-limited 
filled aperture with a diameter of D = 12 m operating at a wavelength of A — 500 nm has a central 
response FWHM of 1.22 X/D = 0.01". In contrast, the final resolution of an image reconstructed from 
a set of interferometer measurements is not uniquely defined. By varying the weighting of the individual 
(it,i>)-coordinates, or even removing some baselines completely, one can put more or less emphasis on certain 
spatial frequencies in the reconstruction, and thus change the effective resolution, even after the data are 
taken. 

For example, using a Gaussian weighting which falls to 0.25 at the longest (12 m) baseline results in 
a FWHM of 7 mas with very low sidelobes. Another option is to use only a narrow ring of long baselines 
which yields a resolution almost twice as high as that of the full telescope. This fact is regularly exploited 
in "aperture masking" observations of bright double stars where collecting area is not an issue, but the 
highest possible resolution is required. Therefore, SIM will yield a resolution for synthesis imaging anywhere 
between 5 and 10 mas, depending on the actual distribution of measurements in the (u,?;)-plane, and the 
weighting algorithm used. 



2.2.2. Field of View 

The field of view (FOV) is limited by the coherence length of the light, I = ct = cj Av = A 2 /AA, where 
r is the coherence time of the light, c is the speed of light, and Av is the bandwidth of the filter used for 
the observations. It can be easily calculated as the product of the spectral and spatial resolutions, which 
can be seen as follows: while the path lengths from the two interferometer elements are assumed to be 
equal in the center of the field, an angular separation O from the center along the baseline vector D results 
in a path difference 8 • D. As long as this difference is small compared to I = A 2 /AA, the fringe contrast 
will not be affected. Thus, < • -£t, or < i? S patiai ■ -Rspcctrai- For example, a 10 nm wide channel at 
A = 600 nm would result in a FOV of 0.6". Any single baseline has a FOV that forms a 0.6" wide, long strip 
on the sky that is positioned perpendicular to the baseline. However, as the baseline rotates, sources along 



2 The exact meaning of "sufficiently" is a non-trivial problem and depends critically on the (unknown) source structure. 
This is the reason why detailed source modeling is needed to predict the imaging performance of any interferometer. 
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the strip more than 0.6" away from the center will appear and disappear in the FOV, thus affecting the 
measurements in an unpredictable way. In order to avoid this confusion, SIM will have a circular aperture 
stop in its optical path. 

There is yet another constraint that sets even tighter limits on the FOV, namely, the need to avoid 
aliasing. Especially for a complicated source structure or crowded field this can be a major limitation on 
any deconvolution algorithm, since the sum of the grating rings of all the sources in the field can — and will 



— mimic fake sources, as will be demonstrated in Section 5A. Unless a priori knowledge about the source 
structure exists, aliasing must be avoided for reliable source reconstruction. This means that in practice 
the SIM FOV is limited to a circle with a diameter of 0.2" , which can be seen as follows. The position 
of the first grating ring of a point source in the field is determined by the spacing of the baselines. In the 
case of 0.5 m increments, the grating ring has a radius of 10/0.5=20 times X/D or 0.2". To make sure that 
no aliasing occurs from any source inside the FOV, it has to be limited to a diameter of 0.2". All of our 
simulations in this paper, however, have been carried out with a FOV of 0.3" in order to demonstrate the 
above effects. 

Figure [j] shows a possible (M,w)-coverage obtained with 170 measurements and the resulting point 
spread function (PSF) or "dirty beam" . Here, and in all the simulations that are described in this paper, 
we have set the longest baseline to 10 m, which corresponds to the requirement of the original SIM design. 
The current design offers a maximum baseline of 12 m. The grating rings in the PSF described above are 
clearly visible. The number of samples in the (u,i>)-plane in Fig. |l| is twice the number of measurements 
made. This is because any astronomical object has a purely real brightness distribution, so the mutual 
coherence function of the incoming wavefront is hcrmitian. Therefore, once the complex visibility of point 
(u,v) has been measured, its conjugate can be entered at point (— u,— v). In addition, different from the 
case of radio interferometry, the total flux of the source is determined during every visibility measurement. 
This gives the value at the origin of the (w,,i>)-plane (the "zero spacing" ) to high accuracy. 

The "noisy" residuals around the central peak of the PSF are due to the incomplete (w,w)-coverage. 
Since, however, the (u,i>)-coordinates of the SIM measurements are defined with extremely high precision, 
the shape of the PSF can be calculated very accurately and without the uncertainties and noise of an actual 
measurement. It is this precise knowledge of the PSF that makes SIM such a powerful imager, since the 
success of deconvolution algorithms like CLEAN or the Maximum Entropy Method (MEM) is limited by 
the knowledge of the PSF. 

The accuracy of the "raw" synthesized image before deconvolution depends in a complicated way on 
several parameters, including the noise properties and flux distribution of the source, the stability of the 
phase measurements, and the actual (u,i;)-coverage. In order to investigate these issues in detail, we have 
developed a software toolkit which we will describe in Section 0. 



2.2.3. Sensitivity 

Jumping ahead, we list here the results of the simulations described in sections ^| and [s] in order 
to provide an estimate of typical integration times needed to perform synthesis imaging of extended 
astronomical objects. 

The surface brightness threshold for making scientifically useful images with SIM in the continuum 
band from 500-800 run is SB C w 0.3 mJy/D" which corresponds to my = 17.7/D". Under the same 
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conditions an emission line such as Ha can be imaged in a narrow (4 nm was assumed in the simulations) 
band with a threshold surface brightness of SBi w 2.5 x 10~ 14 ergs/cm 2 /sec/n". A full synthesis with SIM 
will take typically 4-6 hours of on-target integration time for sources at these brightness levels. During 
this time, 150 - 200 different baselines and baseline orientations will be observed. Observations can of 
course be repeated to increase the sensitivity. 

These sensitivities can be achieved even in the presence of point sources up to V=12 in the held of 
view using the nulling imaging mode discussed in the next section. 

2.3. Imaging in nulling mode 

Because of its unprecedented baseline stability, SIM has great potential for interferometric nulling of 
point sources. One of the SIM beam combiners will introduce an achromatic 180° phase shift in one of 
the inteferometer arms by polarization inversion, thus eliminating the light from a point source that is 
located at the phase center. This capability is a programmatic technology demonstration requirement. 
The experience gained is crucial for the design of future space missions that are aimed at the detection 
of extra-solar planets, like e.g. the Terrestrial Planet Finder (TPF). The requirement for SIM is to reach 
a nulling efficiency of 10~ 4 . Many interesting science programs intend to use this capability to determine 
the spatial extent of objects like stars or supernova shells by measuring the light leakage around the phase 
center. 

In this paper, we will describe another mode of observation, namely synthesis imaging in nulling mode. 
In Section 5, we will discuss the specific features of this technique for the case of exo-zodiacal dust disks. 

2.4. Spectral synthesis imaging 

Since SIM is designed to perform astrometric measurements in a dispersed fringe mode, it has a built-in 
spectroscopic capability. That is to say, the amplitude and phase of an object's wavefront are monitored 
simultaneously at a number of different wavelengths. Since the same physical baseline translates into 
different (w, ^-coordinates at different wavelengths, this information can, in principle, be used to increase 
the (u,u)-coverage. This method is usually referred to as spectral synthesis imaging. 

Spectral synthesis has one major caveat: the implicit assumption in the reconstruction of the source is 
that its structure does not change with wavelength. In the optical regime this assumption is far from safe. 
In fact, it is often plain wrong, especially when strong emission lines are present in the bandpass considered. 
That does not mean that one cannot use spectral synthesis with SIM at all. There are classes of objects 
for which the above assumption is acceptable. For example, in the case of zodiacal disks around solar-type 
stars, it may be reasonable to assume that the structure of the source does not vary strongly over the SIM 
passband, since it is observed only in reflected starlight. We will present the results of spectral synthesis 
simulations for exo-zodiacal disks in Section q. 
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3. SIMsim — a software toolkit 

In order to simulate the imaging capabilities of SIM, we have developed SIMsim, a software package 
written in the Interactive Data Language (IDL) environment. The motivation was to obtain realistic 
performance estimates for a variety of target sources, and to develop optimal observing strategies. The 
principal steps carried out during a simulation are illustrated in Fig. |^. 

In brief, a model source distribution is put on a finely-sampled grid (typically 256 x 256 pixel). For 
simplicity, the field center is assumed to be the center of rotation of the spacecraft during the synthesis 
imaging, and no effort has been made so far to account for possible pointing instabilities. The mutual 
coherence function of the source wavefront is calculated via a Fast Fourier Transform (FFT), and sampled 
on the (u,w)-coordinates observed by SIM. Noise is added, and the actual phase and amplitude retrieval 
process is simulated. The details of this algorithm are described in Appendix B. So far, we have incorporated 
photon noise according to Poisson statistics, and random phase calibration errors with an amplitude of up 
to 0.5°. A final Fourier inversion of this sampled and corrupted data then gives the "dirty" image. The 
PSF is calculated by transforming a mutual coherence function of unit amplitude and zero phase at each 
(u,t;)-sample point of the observation. Deconvolution algorithms like CLEAN or the Maximum Entropy 
Method (MEM) can then be applied to minimize the effects of the incomplete (u,w)-coverage. 

The nulling capability of SIM is modeled as follows. Before the noise is added to the source wavefront, 
its amplitude A (at each (w,z;)-coordinate) is reduced by most of the amplitude A c of the source at the 
phase center u = v = 0: 

Anulled = A- A c ■ (1 - e) (1) 

where e denotes the accuracy to which nulling can be achieved. For SIM, e is expected to be around 1CU 4 . It 
is limited by stray light due to imperfections in the optical elements and by the stability of the pathlcngths. 
The phase of the wavefront is set to zero, resulting in an intrinsically symmetric source reconstruction. The 
reason for the loss of phase information is that the delay cannot be modulated without losing the null: once 
the light is minimized, the delay lines have to be kept as stable as possible. The reduced amplitude of the 
signal results in lower photon counts on the detector which, in turn, reduces the poisson noise. 

Equation [l] ensures that the total flux in the reconstructed image is equal to the flux in the whole 
field-of-view except for the central pixel which, in the simulations, contains the flux of the nulled star. 
This is, of course, known for the models. However, in a real observation, such a priori information on the 
exact value of A c does not exist. It is only possible to measure the total flux - including the central source 
- by doing a standard observation in the direct beam combiner. An iterative approach will be necessary, 
which adds an appropriate constant to each amplitude measurement; this constant is an estimate of 
the total flux in the field, minus the flux from the central source. If this value is chosen correctly, the 
reconstructed image will be zero at the phase center. 

Since the model sources are intrinsically constructed on a regular grid, as are the (u,w)-coordinatcs of 
the "measurements" , there is no need for gridding procedures in our code. However, real sources and real 
observations will not be so conveniently organized, so that in order to fully model the SIM observations, 
weighting and gridding algorithms will have to be implemented. This work is currently underway. 

Figure [| illustrates the Graphical User Interface (GUI) of SIMsim. The top panel allows the user to 
select between a number of model sources and to choose various observational parameters of SIM such as 
(it, v)-coverage, spectral resolution, etc. The intermediate steps of the simulation are displayed in the four 
graphics windows. All results can be saved in postscript format for printing or in FITS format for further 
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work with IDL or any other data analysis environment. In the following sections we use SIMsim to model 
the capabilities of SIM to image two interesting classes of astronomical objects. 



4. Black holes in galactic nuclei 

Black holes have been detected convincingly in at least a dozen galaxies, both active and quiescent. 



The techniques for such detections are multifold (see Kormcndy fc Richstone 1995 for a review), but have 



concentrated on the dynamics of the gas and stars in the neighborhood of the black hole. 

Stellar kinematical studies are often hampered by uncertainties regarding the stellar velocity dispersion 
anisotropy, and the need for high spatial resolution in connection with intrinsically low surface brightness 
makes such measurements very challenging even with HST. However, some galaxies have circumnuclear 
disks of ionized gas which provide a more favorable tool for black hole detection. The emission lines from 
the ionized gas are bright and have large equivalent width, and the dynamical state of the gas is relatively 
simple; its rotation curve directly yields the black hole mass. 

As a starting point for our discussion, we have chosen M 87 as the best known example for such a 
line-emitting gaseous disk. The next section briefly describes our input model for M 87. 



4.1. An example: M 87 



M 87, a dominant giant elliptical galaxy in the Virgo cluster, is the prototype for galaxies that harbour 
a central black hole, the presence of which is revealed by a line-emitting gaseous disk. The nucleus of M 87 



has been imaged with HST in both continuum (Lauer et al. 1992) and line emission (Ford et al. 1994 



and the ionized gas has been studied spectroscopically in great detail (Harms et al. 1994 



Macchctto et al 



1997, Vlarconi et al. 1997) We thus have a fairly clear picture of the properties of the nuclear gas disk. 
Table [if summarizes the parameters we used to construct a realistic model flux distribution for the gaseous 
disk around the central black hole in M 87. The Ha surface brightness over the central 0.25" of our model 
is 1 • 10~ 12 ergs/cm 2 / An", as derived from the observed [SII] flux ( Harms et al. 1994 ) and an assumed line 
ratio of [SII]/(Ha+[NII]) = 4. The gas is assumed to show Keplerian rotation according to 



( 



100 km/s 



f = 0.08869 



Mpc M. 



10 6 M G 



(2) 



where M, denotes the black hole mass, r the angular distance from the center, and d the distance to the 
galaxy. The numbers listed in Table |l| then result in the rotation curve shown in Fig. |I| 

The Ha emission from the rotating disk would appear in a range of spectral channels, according 
to the Doppler shift of the gas velocity projected along the line of sight. For the above parameters, a 
(hypothetical) diffraction-free lOm-telescope in space, equipped with an integral field spectrometer with 
4 nm wide channels, would see the source as in Fig. ||. Here, we have convolved the "true" flux distribution 
of the M 87 model with the "clean beam" , a Gaussian fitted to the central part of the PSF. While the 
central, unresolved continuum source appears in all channels, the appearance of the Ha emission depends 
on parameters like channel width, black hole mass, inclination of the disk, etc. 



To demonstrate the power of our analysis tools, we present in the next section what can be expected 
from SIM observations of M 87 and similar sources. 
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4.2. Results of the SIMsiMulations 

In Fig. we present a few examples of our simulations that are intended to give an impression of what 
imaging with SIM can do in this context. 

The first row of the left panel again shows those channels of our M 87 model containing Ha emission. 
The channels again are 4 nm wide and centered at 554, 558, and 662 nm, respectively^ The second 
row shows the result of the direct Fourier inversion of the complex visibilities measured by SIM. The 
(it, i>)-coverage was the same as in Fig. 1, and the integration time was 100s per measurement, for a total 
on-source integration time of about 5 hours. The nominal SIM throughput of 0.3 was assumed. We added 
photon noise to the measurements according to Poisson statistics. The third row shows the results after 200 
iterations with the MEM algorithm. The emission structure of the high velocity gas is well recovered. The 
same is true for the Ha intensity profile, as shown in Fig. |?]. 

The spatial extent of the Ha emission in each channel provides valuable information on the velocity 
gradient, and thus the black hole mass. This is demonstrated in Fig. |], where we have run the same 
simulation as described above, but for various black hole masses M # . We show only the blueshifted wing of 
the Ha emission for values of M. = 2, 3, and 4xl0 9 M Q . The SIM observations are very sensitive to the 
value of M. , a direct consequence of the linear relationship between black hole mass and radial distance 
in a Keplerian rotation curve, cf. eq. |^. However, it should be stressed that in practice fitting procedures 
are required to solve for all the free model parameters including black hole mass, disc inclination, and the 
actual distribution of disk surface brightness. 

We emphasize that no attempt has been made to optimize the deconvolution process. For example, 
the stability and convergence of MEM can be greatly improved by using "a priori" knowledge about the 
source structure. However, in order to better be able to compare the results for different source models, we 
have — for all the simulations in this paper — started the MEM deconvolution with an empty array, and 
stopped after 200 iterations. For the real SIM observations, one would make every effort to optimize the 
deconvolution strategy for each individual dataset. This is, in fact, another area where future simulation 
work with SIMsim will be useful. However, it is our experience that if no trace of the source structure is 
evident in the direct Fourier inversion, no improvement can be expected from deconvolution, regardless of 
the algorithm employed. Thus, the quality of the data can be most directly judged from the direct Fourier 
inversion itself, before any restaoration. 



4.3. Is there more than M 87? 

An important goal of SIMsiM is to define requirements on various aspects of the SIM spacecraft. As 
an example, what is the minimum spectral resolution of the SIM beam combiner that allows observations 
such as the one discussed here? Certainly, the channel width must not be smaller than the wavelength shift 
of the emission line under observation at a distance of at least one resolution element from the center. This, 
of course, depends on the black hole mass and the distance to the host galaxy. In Fig. |9|, we have plotted 
the limiting black hole mass (i.e. the minimum mass that causes a shift of the Ha emission of at least one 
spectral resolution element at an angular separation of two spatial resolution elements (0.02") from the 



3 We have masked the central continuum source in these plots for better display only; the observation was not simulated in 
nulling mode. 
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center) as a function of the distance to the host galaxy for a number of possible SIM spectral resolutions 
and disk inclinations. For example, with 1 nm wide spectral channels (dotted lines), SIM can potentially 
detect 1O 8 M black holes up to a distance of w 80 Mpc, if the disk has an inclination i < 20°. Whether 
the Ha-flux in cases different from M 87 is sufficient to acquire a high signal-to-noise ratio in a reasonable 
amount of time is more uncertain. 



5. Exo-zodiacal disks 



Circumstellar dust disks are commonly thought to provide the material for planet formation. The 
structure of such disks is therefore important for our understanding of the processes that lead to the 
formation of planetary systems. Many searches have been started for such disks, most of which are based 
on coronographic observations and concentrate on stars with a high IR-excess as measured by the IRAS 
satellite (e.g. [Smith, Fountain, fc Terrile 1992 ). Until most recently, however, only one unambiguous 
detection of a circumstellar dust disk outside the solar system had been reported, namely the (3 Pic system. 



From modeling the (3 Pic disk, Kalas fc Jewitt (1996) found that it is the uncommonly large scattering 
cross-section of the dust particles that makes f3 Pic a special case. If there were other nearby systems like 



(3 Pic, they would have been detected regardless of the disk inclination. Kalas & Jewitt (1996) conclude 
that at least an order of magnitude improvement in the suppression of the stellar light is needed to detect 
more "normal" circumstellar disks, i.e. with scattering cross section less than a tenth of that in (3 Pic. 

SIM, with its interferometric nulling capability, differs from groundbased coronagraphic observations 
in a number of ways: 



1. The starlight suppression in SIM is not limited by atmospheric seeing or the efficiency of an adaptive 
optics system, but only by the quality of the optical elements and the pathlcngth stability. 

2. The high spatial resolution of SIM allows the detection of disks much smaller than that in (3 Pic, 
and SIM can therefore image such systems to greater distances d. In fact, SIM's ability to detect 
exozodiacal disks increases with d, as long as the disk fills the 0.2" FOV of SIM. This somewhat 
surprising behavior is due to the fact that the surface brightness of the disk varies only slightly with 
d, so its total signal on the SIM detector will stay more or less constant. The signal from the star, on 
the other hand, decreases as 1/d 2 , so that the noise contribution from the nulling residuals of 10~ 4 
become less important. In addition, the further away the star is, the more it resembles a true point 
source. This will reduce light leakage around the null, which is also a source of photon noise even for 
a perfect nulling beam combiner. 

3. Since all phase information in the nulling process is lost, the source after reconstruction always 
appears symmetric. That means that eventual asymmetries in the disk can not be recovered. 

4. Aliasing due to the incomplete sampling of the mutual coherence function can sometimes mimic 
companion stars or planets. 



In what follows, we investigate the SIM potential for the detection of zodiacal disks and specifically address 
the points listed above. 
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5.1. The model 

We use a simple model of a circumstellar dust disk, based on the zodiacal emission of the solar system 



(Reach et al. 1996, Traub 1998). The particle number density of the solar system zodiacal dust cloud can 
be well modeled by the expression 

p = po ■ r- a ■ e- fHz '^ (3) 

where r is the in-plane radial distance from the center, z is the modulus of the perpendicular-to-plane 
distance (in AU), a = 1.4, (3 = 3.26, and 7 = 1.02. We treat the scattering efficiency of the dust particles 
as a flux scaling parameter, and we assume isotropic scattering. In what follows, we investigate the two 
limiting cases of a face-on (?=90°) and edge-on (i=0°) orientation of the disk. 

The radial surface brightness distribution of the face-on disk can be derived from integrating the 
volume density. For 7 = 1 it is given by the expression 

t7 ( <=90o)=C-criT- o ' +1 . (4) 

Here, C is a scaling factor that characterizes the scattering efficiency of the dust particles in an exo-solar 
system relative to the solar system. The normalization constant a\ is the surface brightness of the face-on 
disk, as seen from outside the solar system, at a distance of 1 AU from the sun. Assuming isotropic 
scattering, this number is (at A=550 nm) 

0-1 = 0.13 MJy/sr. (5) 

An approximate ratio of the surface densities of the edge-on and face-on directions at a given value for r 
can be obtained by ratioing the characteristic scale heights of the z-direction and along a tangent to r in 

the plane: 

e/f = 0>-h . Ve 2 /" - 1 « 5. (6) 

An approximation to the edge-on surface brightness distribution is therefore 

cr (i=0 c ) =C-e//-a 1 T- a + 1 -e^ £ . (7) 

where 7 « 1 is used. Figure |h] gives a visual impression of the model in both the face-on and edge-on case. 
The spectrum of the dust in the wavelength range discussed here (500 - 1000 nm) is assumed to be that of 
a black body with the temperature of the central star. 



5.2. j3 Pic — the showcase 

In this section, we investigate the question whether SIM is capable of finding other systems like (3 Pic, 
should they exist. In this context, we define "/3 Pic-like" as having about 1000 times the dust content and 
surface brightness of the solar system, i.e. C=1000, without any detailed modeling of the specific properties 
of /3 Pic itself. 

As mentioned above, the SIM performance in terms of imaging of circumstellar disks is very limited 
for nearby systems, but improves with increasing distance, as long as the disk fills the FOV. For example, 
the (3 Pic disk itself is impossible to image with SIM, because even with 10 4 nulling, the photon noise of 
the remaining light from [3 Pic would make the detection of the faint disk emission impossible. Figure [ll] 
demonstrates the increasing sensitivity of SIM up to distances of 1 kpc. 
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Again, we note that no attempt has been made to optimize the restoration procedure in order to avoid 
additional degrees of freedom. As can be seen in Figure [n], SIM can easily detect (3 Pic-like systems up to a 
distance of about 1 kpc, beyond which our model disk no longer fills the SIM FOV and its total flux drops 
rapidly. This distance limit depends linearly on the extent of the disk, for which we have assumed a value 
of 100 AU. This is a somewhat arbitrary number, chosen such that the angular diameter at a distance of 
1 kpc is 0.2", equal to the SIM FOV. The central star is assumed to be of solar type (My— 5). 

5.3. Can SIM find other solar systems? 

The above results change dramatically if systems with lower dust contents than (3 Pic are considered. 
The dust content of our model disks is parametrized by the constant C in eq. ^. Figure [l^ summarizes our 
results on studying various values for C. Disks with more than a tenth of the dust content of (3 Pic can be 
detected with SIM. However, it is evident that zodiacal disks similar to that of the solar system will not be 
found with SIM, regardless of their distance or inclination. 

5.4. What are the clumps? 

It is interesting that in all of our reconstructions, the disk structure after deconvolution appears much 
more clumpy than the smooth input model. In fact, when looking at the results, one might in many 
cases think to have found a number of planets. These "clumps" are artefacts of the restoration process 
on our incomplete (u,i>)-data. While the sidelobe pattern itself is known to high accuracy, the poisson 
noise associated with the signal does not have the shape of the PSF. Spikes in the noise distribution create 
residuals that are mistaken as real source structure by any deconvolution algorithm. The fact that the 
clumps are indeed randomly distributed can be seen when comparing three independent simulations of the 
face-on zodiacal disk of a (3 Pic-like system at 1 kpc and a nulling efficiency of e = 10~ 4 in Figures |IT| , |l2| , 
and |l^, which only differ in their random photon noise. A repetition of the observation would allow 
identification of such spurious objects. 

5.5. How critical is the Nulling? 

In order to investigate the effects of a degradation of the SIM nulling capability, we have run a series 
of simulations with different values for e, the nulling efficiency in eq. [I]. Figure |l^ demonstrates that for a 
(3 Pic-like system at 1 kpc, the most favorable case in all of the simulations we have run so far, SIM could 
detect the disk even if it had only e = 10 -3 . However, anything worse than that would not suppress the 
starlight sufficiently even at such high distances. We conclude that in order to be a powerful search engine 
for exo-zodiacal disks, SIM needs to achieve at least e = 10~ 4 , and anything better would improve the 
performance considerably. 

5.6. Planetary Gaps 

A massive planet — if present — is believed to sweep up any inward-spiraling debris and thus cause 
a dust-free region in the circumstellar disk inward of the planet's orbit. Such an "inner gap" is believed 
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to exist in the (3 Pic-system, as inferred from its spectral energy distribution (Lagage & Pantin 1994), but 
the gap has never been directly imaged owing to the large radius which is masked out by coronographic 
instruments. 

With its high spatial resolution and sensitive nulling, SIM is able to directly image possible inner clear 
regions in exo-zodiacal systems. This is demonstrated in Fig.[l4] where we have investigated the case of a 
j3 Pic-like disk at a distance of 500 pc with and without a 20 AU inner clear region. Clearly, the signature 
of the disk depends on its size and the distance to the system, but a 20 AU radius gap would be detectable 
in systems as distant as 1 kpc. Here, we just want to illustrate the intriguing prospect of directly imaging 
such gaps (and later, possibly the planets that cause them). 



6. Summary 

We have built an interactive software package that enables us to perform numerical simulations of 
imaging observations with SIM. The results obtained so far can be summarized as follows: 

1. SIM, with a spectral resolution of 1 nm, can resolve the velocity structure of line-emitting gaseous 
disks around black holes in the centers of nearby galaxies down to black hole masses of 10 8 M Q at 
distances up to 100 Mpc. With this resolution, SIM is in principle capable of determining the mass of 
the black hole to 10% accuracy. The necessary integration time depends on the Ha-flux from the disk; 
for cases like M 87, about 5 hours is sufficient. A coarser spectral resolution of SIM would constrain 
the accessible parameter space for quantities such as distance, black hole mass, or disk inclination. 

2. Given a nulling efficiency of 10~ 4 , SIM can image zodiacal disks with dust contents similar to that of 
the j3 Pic-system with high signal-to-noise ratio in about 5 hours of integration time. The limiting 
distance depends only on the angular size of the disk, which should be no smaller than 0.1". 

3. Inner clear regions in the disk can be directly detected and imaged with SIM for systems at a distance 
of up to 1 kpc provided that the dust-free region has a radius of at least 20 AU. 

We are grateful to R. van der Marel for his advice on the black hole simulations, and to W. Traub for 
his expertise on the zodiacal disk model. We also would like to thank S. Unwin, M. Colavita and M. Shao 
for their interest in and continued support of this project. The work described in this paper is carried out 
at STScI with financial support from the SIM project at the Jet Propulsion Laboratory (JPL). 



A. The Fourier optics of SIM 

The idea that an interferometer could be constructed to measure the angular diameters of stars is 
usually credited to Fizeau (1868), although no details were provided in that particular reference. Twenty 
years later, Michelson (1890) took it up and, another 30 years later with Pease, constructed the now-famous 
stellar interferometer on the 100- inch telescope ( Michelson fc Pease 1921 ). Figure [l5| shows a sketch of the 
original drawing from that paper, to which we have added symbols indicating the spacing L of the outrigger 
mirrors and the spacing D of the mirrors at the telescope aperture which direct the light to the primary 
mirror. Variations of this sketch appear in many introductory textbooks on physics and astronomy. 



- 13 - 



The image of a star in the focal plane of the stellar interferometer is crossed by interference fringes, 
as shown in cross section in Fig.|l6|, where some additional parameters are defined: the fringe amplitude 
A, the fringe period A9 and the fringe phase <p. In this picture the average signal from the star has been 
normalized to 1. The overall size of the Airy pattern is set by the aperture formed by the cascade of small 
reflecting mirrors which direct the light to the telescope primary mirror. 

As is well known from the general theory of wave diffraction, this Fraunhofer diffraction pattern in the 
focal plane of the telescope is the square modulus of the Fourier Transform of the distribution of illumination 
(in amplitude and in relative phase) over the aperture. The expression for the stellar interferometer 
can be derived in a straightforward way, at least in one dimension, by starting from the result for the 
single-slit diffraction problem. We then use standard theorems from the theory of Fourier Transforms, in 
particular the shift theorem and the convolution theorem, and apply them to a plane wavefront over the 
entrance aperture (pupil) . In this way it is straightforward to show that the final Fraunhofer pattern for an 
unresolved star is the square modulus of 

M{u) = F(u-u )xG(u-ui) (Al) 
= 2 sinc(7ra(u — Un)) cos(ir(Du — Lug)) 

with sinc(a;) = sin(a;)/2;, and u = sinifi, <\> being an angle as defined in Fig.|l6|. This is the product of 
the pattern F(u) — sinc(7ra(u — uo)) for a single slit of width a offset to uq = sin</>o/A, and the pattern 
G(u) = 2 cos(ttD(u — Luq/D)) for a very thin double slit offset to u\ = Luq/D, where 4>o is the angular 
distance of the star from the direction normal to the baseline L. We note that: 



• As the source position cf>o changes, the fringe phase moves at an amplified rate Luq/D. The positional 
accuracy is therefore magnified by L/D. 

• The fringe period remains set by the separation D of the secondary apertures no matter what is the 
separation L of the primary apertures on the outriggers. 



SIM can be regarded as a Michelson stellar interferometer with two specific features: 

• The separation D of the secondary apertures is reduced to zero^j. The fringe period therefore becomes 
infinitely large, and the entire image of the star is either light or dark depending on the fringe phase. 

• A delay A = L sin <f>Q + S is added to one side of the interferometer light path (e.g. between Mi and 
M 2 in order to shift the double-slit pattern from u\ back to within a small offset 5 from the position 
of the star u a . This is necessary because of the finite bandwidth of the signal, which results in a 
limited coherence length, as discussed in Section 2.2.2] . We will discuss the bandwidth effects on the 



fringe contrast later in this section. 

The Fraunhofer pattern for a monochromatic point source in our one-dimensional model of SIM is therefore 
the square modulus of: 

M(u) = 2F{u-u ) cos(vr(L sin </>o - A)/A), (A2) 

which is 

P(u) = \M(u)\ 2 = P (u) (1 + cos(27ntyA)) (A3) 



4 Practical methods of achieving this beam combination usually introduce an additional fringe phase of typically fa it/A which 
we ignore here. 
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wherc Po(u) = 2 \F(u — uo)| 2 - Viewed as a function of the delay offset S, this is a co-sinusoidal pattern with 
an average value set by the intensity Pq of the star^j, and varying from zero to 2P$. The fringe amplitude 
here is 1, appropriate for a point source. 

Consider now an extended source in the sky, which we integrate in a direction perpendicular to the 
interferometer baseline in order to produce a 1-D image B{(f) where <f> is an angle measured in the plane 
of the baseline (Fig. [Tt]). Let us further assume that the position cf>o to which we have initially offset the 
fringe pattern from our calibration grid star (by adding a delay A) is somewhere in the middle of B((f>). 
Each part d<p of the spatially-extended source then contributes a fringe pattern: 

dP{<p) =S(0)#(1 + Sft e -2™(£sm0-A)/A)). ( A4 ) 

where again A = Lsin^n + 5 and 5ft denotes the real part of the exponential. Let us further simplify the 
notation by assuming that the FOV is confined to a small region on the sky, so that sin0 w </>, and measure 



L in wavelengths, I = L/X. The total fringe pattern obtained by integrating equation A4 is then: 

r + co 

P = B Q A(t) + ^ B^e^^-^^'^d^. (A5) 



We can write the integral in this equation as: 

B A(f> X SR e 2«(i0o+5/A) x ( A g) 

where we recognize T(l) as the Fourier transform of the normalized B((f>): 

r(z)= 00 tT a / -• (A7) 

B A(p 

Since B((f>) will in general not be symmetric about </>o, T(l) will be complex, so we can write 
T(l) = A{l)e 2m ^ l \ leading to: 

P = P a {l + A(l) cos(2Tr(l(j) - cj)(l) + 5/X))}. (A8) 

There are 3 unknowns here: the average power Pq, the fringe amplitude A(l), and the fringe phase <p{l) 
(measured here in turns). They can all be determined by measuring P as a function of delay 5. In practice, 
this can be achieved by continously scanning 8 over ±A, and recording P over a number of time intervals. 
This is the basis for our simulation of the measurement process and the introduction of photon noise, as 
will be described in appendix B. 

Let us now consider the effects of the finite bandwidth of the interfering light. In order to keep the 
algebra manageable, we consider only the case of a plane wave from a point source. The fringe pattern we 



have derived in equation A3 will be different for different values of the wavelength A. The total response 



is the sum of the fringe patterns for all wavelengths over the spectral range of the passband. Suppose 



the passband has a spectral shape S(v). Then, in analogy to equation A4, that part of the spectrum at 
frequency v contributes a fringe: 

dP(v) = S{u)dv{l + SR e 2^ l (isi„0 o -A)/A^ ( A9 ) 



5 The residual factor of 2 present in equation A2 is ignored further here. It comes in because we have doubled the collecting 
area by using two outrigger mirrors, and it will be subsumed into Pq by the overall interferometer calibration procedure. 
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where Pq = S(v)dv. We define a delay time r — 5/c where c is the speed of light. Then we have: 

dP{u) = S(v)dv{l + Ue- 2 ™ T }. (A10) 
The total fringe pattern is obtained by integrating: 

/+oo 
S{v)e~^ lvr dv. (All) 
-oo 

We recognize the integral as the Fourier transform s(t) of the filter function S(v). As an example, let's take 
a rectangular filter of height So and width Av centered at vq. The transform of this function is related by 
the shift theorem to the transform of a function C(v) having the same shape but centered at the origin. The 
transform of C(y) is c(r) = So /S.v sinc(7rrZ\^), so the transform of S(v) = C(v — vq) is s(r) = e 27rwor c(r), 
and the fringe pattern becomes: 

P{u) = 5 A^{1 + s\\\c(tttAv) cos{2ttv t)}, (A12) 

so that if we scan 8 = ct and record the intensity of the stellar image in the focal plane of SIM we will 
obtain the pattern shown in Figure [181 



B. Simulation of the measurement and its noise 



In order to provide an estimate of SIM sensitivity and image quality various noise components 
associated with the measurements need to be modeled as realistically as possible. The systematic 
(calibration) errors of the spacecraft are hard to predict at this time, because the design has not been 
finalized and a detailed hardware model is not yet available. As of now, we can only make assumptions on 
parameters such as fringe phase stability. On the other hand, the statistical errors due to photon statistics 
are straightforward to model. They are implemented in SIMsim as follows. 

The fundamental shape of the signal as a function of path difference (delay) between the two interfering 



beams has been derived in equation A8. Here, the term "signal" means the total number of photons falling 
on the detector over the integration time. If the delay <5 is continously modulated over ±A, then one can 
integrate the detector signal over four time "bins" , during each of which the delay moves over a quarter 



wave. This method has first been described by Wyant 1975 and was suggested for use in astronomical 
interferometers by Shao & Staelin (1977). It is usually referred to as the "quadrature method" and has 



been used in many ground-based interferometers. The usual implementation uses a triangular modulation 
of the delay with piezo-electric elements (e.g. [3hao et al. 1988 ). 



In ground-based interferometers the delay is scanned over the range of one wavelength in about 10 ms 
owing to the short coherence time of the atmospheric seeing, which otherwise moves the fringe phase 4> in 
a random way. In principle, there is no need for high modulation frequencies in space, so that instead of 
a continous modulation, one could envisage a "stop-and-stare" mode for the delay 5. However, since the 
continous modulation technique is well established and simple to handle mathematically, we have used it to 



simulate the measurements in SIMsim. Following [Wyant 1975| , the signal in the four bins B1...B4 which, 
for mathematical simplicity have the limits [—A/8, A/8], [A/8, 3 A/8], [3 A/8, 5 A/8], and [5A/8,7A/8] can 



calculated by integrating equation A8 over these intervals, yielding 



Bl u>v = P /4+ 1/(ttv / 2) A{u,v) sm<j){u,v); 
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B'S 



B2 



B4, 



P /4+ 1/(ttV2)A(u,v) coscj)(u,v); 
Po/4- l/(wV2) A(u,v) sin0(w, v);and 
P /4-l/(irV2)A(u,v) cos<p(u,v). 



Po here is the total signal (in photons) in the FOV over the spectral passband during the integration time. 
In the case of a simulation, the values for the fringe amplitude and phase, A(u,v) and <p(u,v), are known 
from the Fourier transform of the model source structure. Therefore, the values Bl — B4 can be calculated, 
and poisson noise can be added to them, before using the "measurements" to recover the source. 

The noisy bin signals are then used to determine the "measured" fringe amplitude and phase according 

to 



These values are re-entered in an empty array (the (u,v)-plane), and the (sparsely-filled) array is 
Fourier back-transformed to yield the reconstructed image of the source. As an option for the user, an 
additional noise component can be added that tries to model the expected phase calibration errors. Our 
simulations have shown, however, that up to the specified maximum error of 0.5°, this has negligible 
effects on the restoration quality. Since the PSF can also be computed, the recovered image can further be 
improved by applying restoration algorithms like CLEAN or MEM. 



A 



noisy 



7T/V2 x ^{Bl - B'S) 2 + (B2 - 54) 2 



noisy 
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Fig. 1. — Left: example of a specific (u,i>)-coverage, obtained from 170 SIM measurements. Right: the 
resulting point-spread-function. 

Fig. 2. — Schematic view of SIMsim capabilities . 

Fig. 3. — Graphical User Interface of SIMsiM. 

Fig. 4. — Model rotation curve for M 87. 

Fig. 5. — Model surface brightness distribution after convolution with the "clean beam" for the Ha line 
emission from a rotating gas disk around the black hole in M 87. The channels are 4 nm wide and centered 
at [-8,-4,0,4,8] nm offset from the line center. 

Fig. 6. — Results of the simulated SIM synthesis imaging observations of the nucleus of M 87. The integration 
time was 100 s per baseline, and we used 170 baselines for a total on-source integration time of 5 hours. 
Shown are contour plots of the model source after convolution with the "clean beam" (top row), the direct 
Fourier inversion (center row), and the deconvolved Fourier inversion (bottom row) after 200 MEM iterations. 
A 0.05" region around the central bright continuum source has been masked after the simulation to better 
display the extended stucture. Contour levels in all plots are 10, 20, 30. . . 100% of the peak surface brightness, 
which is given below each plot in units of 10 -13 ergs/s/cm 2 /A/n" 

Fig. 7. — Profile plots for the model source (top), the direct Fourier inversion (middle), and the restored 
image after 200 MEM iterations (bottom). 

Fig. 8. — Structure of the blueshifted wing of the Ha emission (top), direct Fourier inversion (center), and 
reconstructed image after 200 MEM iterations (bottom) for black hole masses of 2, 3, 4-10 9 M Q . Contour 
levels are as in Fig. |6| . 

Fig. 9. — Limiting black hole mass for the detection of a rotating gas disk around a black hole with SIM. 
The Doppler shift of the line has to be at least equal to the SIM channel width. This limit is shown as a 
function of black hole mass, disk inclination, and SIM spectral resolution. The solid lines refer to 4 nm wide 
SIM channels, dotted lines are for 1 nm channels, and dashed lines for 16 nm channel width. Four lines per 
set are shown, referring to disk inclinations of i=0° (edge-on), 20°, 40°, and 60°, the last being the highest 
curve of each set. 

Fig. 10. — Surface brightness distribution of our zodiacal disk model, in face-on (left, i = 90°) and edge-on 
orientation (right, i = 0°). The grey scale is linear and identical in both plots to emphasize the higher 
surface brightness in the edge-on case. The central star has been removed for better contrast in the display. 

Fig. 11. — Results of the SIMsiM reconstruction for a (3 Pic-like zodiacal disk oriented at an inclination 
of 0° (edge-on, left panel) and 90° (face-on, right panel). The columns in each panel are for distances of 
100 pc (left column), 500 pc (center), and 1 kpc (right). The rows in each panel are for the input model after 
convolution with the "clean beam" (top row), the direct Fourier inversion (center row) and reconstructed 
image after 200 MEM iterations (bottom row). A 0.05" diameter region around the central star has been 
masked out after the simulation to better show the extended, low-level structure. All panels comprise a 0.3" 
field of view. Contour levels are 10,20,30. . .100% of the maximum (unmasked) surface brightness in each 
image (given below each panel in units of 10~ 15 ergs/s/cm 2 /A/D"). The assumed SIM nulling efficiency is 
e = 10~ 4 . Data were collected in five spectral channels (100 nm width between A = 500 nm and A = 1000 nm) 



and used for spectral synthesis imaging. 



Fig. 12. — Results of the SIMsim reconstruction for models with various dust content in the case of an 
edge-on (left) and face-on disk (right). The columns in each panel denote (from left to right): C=1000, 100, 
10 times the solar dust content. Field size and contour levels are as in Fig. 

Fig. 13. — SIMsim reconstructions of a /3 Pic- like system (C=1000) at a distance of 1 kpc. Shown are the 
results for various values of the SIM nulling efficiency e. Again, the edge-on case is shown in the left panel and 
the face-on case on the right. The columns in each panel denote (from left to right): e = 10 -4 , 10~ 3 , 10 -2 . 
Field size and contour levels are as in Fig. [ll]. 

Fig. 14. — SIMsim reconstruction of a /3 Pic-like zodiacal disk (C=1000) at d —500 pc without (top) and 
with (bottom) a 20 AU radius inner clear region. The central star has been masked after the simulation to 
better show the faint disk emission. 

Fig. 15. — Sketch of Michelson's stellar interferometer on the 100-inch telescope. 

Fig. 16. — Focal plane profile of a star in Michelson's stellar interferometer. 

Fig. 17. — Geometry for collapsing a 2-D source onto the fringe normal direction. 



Fig. 18.— 



SIM fringe pattern for a rectangular filter function. 
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Table 1. Input parameters for the M 87 model 



Property Ref. 

stellar background /ii = 15.0 mag/d" 1 

central point source I = 16.1 1 
H Q (+ [Nil]) disk flux 1 ■ 10~ 12 ergs/cm 2 /s/D" 2 

Radial intensity profile Z(r) = A ■ e r / 01 

black hole mass 3 • 10 9 M Q 3 

distance 16 Mpc 

inclination 45° 

position angle 0° 

grid resolution 0.005" 



Note. — References: 1) Laucr et al. 1992| 2) Harms e\ 
al. 1994 assuming a line ratio of [SII]/(Ha+[NII]) = 4; 3) 
|Marconi et al. 1997 . 
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